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I. 


INTRODUCTION 


The Supersonic Gas-Liquid Cleaning System Research Project consisted mainly of a 
feasibility study - including theoretical and engineering analysis - of a proof-of-concept 
prototype of this particular cleaning system developed by NASA-KSC. The cleaning 
system utilizes gas-liquid supersonic nozzles to generate high impingement velocities at 
the surface of the device to be cleaned. The cleaning fluid being accelerated to these high 
velocities may consist of any solvent or liquid, including water. Compressed air or any 
inert gas is used to provide the conveying medium for the liquid, as well as substantially 
reduce the total amount of liquid needed to perform adequate surface cleaning and 
cleanliness verification. This type of aqueous cleaning system is considered to be an 
excellent way of conducting cleaning and cleanliness verification operations as 
replacements for the use of CFC 1 13 which must be discontinued by 1995. 

To u tiliz e this particular cleaning system in various cleaning applications for both the 
Space Program and the commercial market, it is essential that the cleaning system, 
especially the supersonic nozzle, be characterized for such applications. This 
characterization consisted of performing theoretical and engineering analysis, identifying 
desirable modifications/extensions to the basic concept, evaluating effects of variations in 
operating parameters, and optimizing hardware design for specific applications. 

Under a partnership arrangement, TRDA, NASA-KSC, Bethune-Cookman College, 
Embry-Riddle Aeronautical University, and Precision Fabricating and Cleaning, working 
as a team, provided the research and technology necessary for the further design and 
development of an operational cleaning system based on the supersonic nozzle concept for 
cleaning and verification of cleanliness. As a commercial partner, PFC plans to utilize this 
research and technology to possibly adapt this cleaning concept for other potential 
markets. 


n. TRDA PARTNERSHIP TASKS 


As per Cooperative Agreement No. NCC 10-00 10, the tasks to be performed by TRDA 
are listed below, along with a brief narrative of activities undertaken by TRDA to satisfy 
our partnership responsibilities. Descriptions provided below are for activities subsequent 
to those documented in the project Performance Report, submitted in January, 1995. 

(1) Provide co-funding (in-kind contributions) of the project. 

TRDA staff spent considerable time and effort managing and coordinating 
project activities. Contributions attributed to the project were in excess of 
$15,400. 



(2) Provide the necessary manpower for proper implementation, 
and business and technical management of the overall project. 

Frank Kinney (TRDA Executive Director) served as Project Director; 
Matthew La Vigne (TRDA Contract Manager) served as Project Manager, 
and Libbie Strickland (TRDA Executive Secretary) provided the required 
secretarial support. The resumes of Kinney and La Vigne were provided 
with the project proposal to NASA. TRDA also utilized qualified 
consultants for technical support on the project. 

(3) Contract with Bethune-Cookman College (B-CC) to provide the 
required technical assistance from qualified engineering students, and 
to assist the TRDA in regard to contracting with E-RAU for fluid 
mechanics experts and managing certain administrative functions, as 
needed. 

The contracting process and roles of each partner were described in the 
Performance Report. 

(4) Serve as the main interface and point of contact for all of the 
project partners. 

TRDA served as the central point of contact for all partners relative to the 
project, while encouraging direct contact and working relationship between 
all the partners. TRDA scheduled and hosted monthly meetings of all 
partners to discuss project activities and review schedules, as well as to 
provide an open forum to consider project obstacles encountered and 
overall direction of the project. 

TRDA also held many informal meetings with individual partners, 
acting as an intermediary on numerous occasions to facilitate working 
relationships between project participants. 

(5) Conduct quarterly progress reviews of the work performed by 
B-CC and E-RAU in conducting the research. 

Although promising only quarterly progress reviews of the work performed 
by project partners, TRDA has took an even more aggressive position in 
overseeing project activities. By contract, TRDA required monthly written 
progress reports from each partner , in addition to a monthly meeting. 

Monthly reports received since the submission of the Performance Report 
from each project partner are included in this report as Attachment “A”. 
These reports give an update on project activities from January through 
March, 1995. 



(6) Track all performance schedules and insure adherence thereto. 

Project schedule information was included in the Performance Report and 
is touched on in the Final Reports contained within this report. 

(7) Provide best efforts in insuring the protection of the partners’ 
proprietary data and intellectual property rights. 

As with all TRDA partnership projects, continued to exercise its best 
efforts to insure the appropriate protection of proprietary data and 
intellectual property rights of the partners. 

(8) Acquire and submit reports and other documentation regarding 
the overall project 

TRDA submitted a Performance Report in January of 1995, and hereby 
submits a Final Summary Report. Activities and findings of the entire 
project are described in great detail in the Final Project Reports from each 
of the three cooperating research partners; Bethune-Cookman College, 
Embry-Riddle Aeronautical University, and Precision Fabricating and 
Cleaning, Inc. These Final Project Reports are included in this report as 
Attachment “C”. 

(9) Prepare and submit the required invoices to NASA-KSC, and 
maintain adequate files and records in connection the cooperative 
agreement and other activities regarding the partnership. 

TRDA submitted the appropriate Federal Cash Transaction Reports to 
NASA-KSC for to cover associated costs of the project. A summary 
outlining the expenditure of funds can be found in Attachment “B” of this 
report. The project expended $ 1 6, 1 0 1 . 54 less than originally budgeted, 
and TRDA requested that NASA-KSC de-obligate these funds from the 
project. 

TRDA maintains a complete file of project related documents, activities, 
etc. relevant to the Supersonic Gas-Liquid Cleaning System. This file is 
open to NASA-KSC inspection at any time. 


ATTACHMENTS 


“A” - Progress Reports 

“B” - Final Project Budget Summary 


“C” - Final Project Reports 



ATTACHMENT “A” 


Supersonic Gas-Liquid Cleaning System Research Project 
Monthly Progress Reports 



SUPERSONIC GAS-LIQUID CLEANING SYSTEM 

MONTHLY PROGRESS REPORT 
January 12, 1995 


This project report is submitted by Dr. Sunil K. David to Bethune Cookman College, 
and describes work performed during the fourth month of the research project (Ref: TRDA 
Grant #411). 

During this month of Dantec 55M01 main unit with cover plate and 55M10 Standard 
Bridge was received at KSC/NASA from the LaRC/NASA. According to the work plan 
mentioned in earlier reports, this instrument is to be set at a convenient place for the velocity 
measurement required to verify the theoretical model. Considering various factors PF&C, 
cocoa was chosen for the site for experiment. Mr Tom Jezowski set all the necessary facility 
for the Nozzle and the control panel in a mobile cleaning lab at PF&C which helped 
considerably in setting the experiments. 

With SGLCS in operation, following task were carried out while: 

1. Anemometer required a probe to measure supersonic velocity of air and air- water mix 
(Gas-Liquid mix) stream coming out of the nozzle. A wedge-shaped hot-film probe (Dantec 
probe 55R31) was selected and procured from Denmark on urgent basis to save the time. 

2. Probe was mounted as per prescribed procedure on a calibrated linear bench. 
Preliminary experiments were performed to put CTA and the Probe in operation. 

3. A procedure was evolved to align the probe, nozzle and the direction of stream in a 
straight line for all the positions of the probe on the linear bench. 

4. Experiments were performed to measure the velocity of the air-stream at various 
distances from the exit-plane of the nozzle for back pressures 200, 300, 400, and 500 PSIG. 
Consistency and the reproducibility of the results were also checked. 

5. Observations under the same conditions will be taken for air-water mixtures to 
investigate the physical effects of the water droplets (or ice dust) in the system.. 

To calibrate the probe either a known flow source is required or a conical pitot will 
be required. 



SUPERSONIC GAS-LIQUID CLEANING SYSTEM PROJECT 

MONTHLY PROGRESS REPORT 
February 10, 1995 


This project report is submitted by Dr. Sunil K. David to Bethune Cookman 
College, and describes work performed during the fifth month of the research project 
(Ref: TRDA Grant #411). 

A meeting was arranged with Dr. James M. Meyers 1 ’ 2 and his colleagues at 
NASA/Langley Research Center, Hampton (Virginia) to discuss methods for particle 
sizing and flow measurement using Laser Doppler Anemometery and the availability 
of the equipment facilities. This visit explored the possibility of using unique facilities 
at Langley and / otherwise to find other economical methods to collect data required 
for the project 3 . Dr. Meyers and his colleague Dr. William M. Humphreys, Jr. have 
kindly agreed to provide their non-intrusive 'Laser Transit Anemometry' facility to 
collect required data on velocity and particle size distribution. 

Experiments were also conducted for several days in the mobile cleaning lab at 
PF&C, Cocoa to measure velocity of the supersonic flow of two phase fluid of the 
nozzle. As mentioned in earlier reports and brief notes, pitot probe is to be used as 
one of the methods to obtain velocity data along the axis of the fluid stream of the 
supersonic nozzle. There are several advantages of using pitot probe for supersonic 
flow measurements 4 . Somehow, a pitot probe for the supersonic flows was not 
available easily, therefore, it was decided to make one. A pitot probe was designed 
at B-CC and tested at PF&C. This probe was used to measure static pressure (after 
the shock) along the axis of the nozzle stream. Readings were taken for the breathing 
air and also for the mixture of breathing air and water at different settings of 
pressures and axial distances. 

Few optical arrangements 5 were also tried this month to get visual images of 
the supersonic flow of the stream of the breathing air and also air-water mixture. 
Finally, after several trials and experimental adjustments clear images of the flows 
were obtained on the screen. Photographs of these images were taken for each 
setting. Equipment from B-CC labs were used for this purpose. 


1. Performance Test of Laser Velocimeter System for the Langley 16-Foot Transonic Tunnel, 

1. James F.Meyers, ec al; NASA technical Paper 2502 /AVSCOM Technical Report 85-B-4, December 1985. 

2. The Elusive Third Component, J.F. Meyers; International Symposium on Laser Anemometry, Ed: A. Dybbs and P Jl. Pfund, Nov. 1985. 

3. Brief note sent to Dr. Meyers regarding the nozzle project and the definition of required data (enclosed) 

4. "Pitot tube as a calibration device for turbulence measurement", Sadek Z. Kassab; Rev. Sci. Instrum. 61(6), 1757-9 (1990) 

5. 'High Speed Wind Tunnel Testing’, Alan Pope and Kennith L. Go in; pp229 ; Robert E. Kriegcr Publi. NY 1978. 
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One of the tasks given to the students working on this project, is to arrange 
and compile the technical literature and subsequently prepare a bibliography of the 
experimental work which will go in the final technical report of this project. This work 
is now in progress and will be finished soon 8 . 

Arrangements are made to conduct phase three experiments 6 7 at NASA/Langley 
research center, Hampton, during Feb. 21-24, 1995. 


6. <"--1— » working on this project could not s tan their work before the Spring session of 1995 due to their earlier commitments for other courses in 

Fall semester schedule of 1994. 


7. ’Schedule of the experimentation, a note submitted in Nov.94. 
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SUPERSONIC GAS-LIQUID CLEANING SYSTEM PROJECT 


MONTHLY PROGRESS REPORT 
March 10, 1995 


This project report is submitted by Dr. Sunil K. David to Bethune Cookman College, 
and describes work performed during the sixth and the final month of the research project 
(Ref: TRDA Grant #411). 

During the last week of February phase-ill experiments were conducted at the 
NASA/Langley research station, to measure of the fluid velocity and the size distribution of 
the particles of air-water mixture in supersonic stream of the cleaning nozzle by a non 
intrusive method. As mentioned in last report, instrument facility of the Laser Transit 
Anemometry (LTA) was made available to us at the LaRC, Hampton, Virginia. The proto- 
type nozzle and the control panel was transported to LaRC before hand. At LaRC all the 
necessary arrangements were made to setup the nozzle and the control panel, and LTA in 
a laboratory which has supply line of high pressure air and Life protection facilities. Tasks 
of optical and Mechanical alignments were also performed using standard techniques. 
Velocity of the particles in two-phase flow of air-water mixture was measured at the various 
axial distances from the exit plane for three upstream pressures (200, 300, and 400 psig) of 
the nozzle. The mass flow rate of water in the air-water stream was also measured at these 
working pressures. Though the ETA method is quite suitable for the absolute measurement 
of the velocity, it is not sensitive enough for determining particle sizes. However, 
experimental results indicated that particle size of the water droplets, on the average, is 
much less than 0.8 micron. 

A detailed description of the above mentioned measurements are presented in the 
final consolidated report which includes tables, spread-sheets of the observations, ambient 
conditions, graphical plots, photographs etc. This report also contains details of other 
methods which were used to obtain experimental data on the performance of the prototype 
nozzle. 

This report will not complete without acknowledging the personal help of Dr. 
William J. Humphreys in setting up the ETA system and collecting data. Thanks are also 
due to Dr. J. Meyers for allowing the use of facilities at Langley and for his expert advice. 




SUPERSONIC GAS-LIQUID CLEANING SYSTEM 

MONTHLY PROGRESS REPORT 


This progress report is submitted by Embry-Riddle Aeronautical Univer- 
sity (ERAU) to Bethune Cookman College (BCC), and describes the work done during 
the fourth month under TRDA Grant #41 1 towards the development of a Supersonic 
Gas-Liquid Cleaning System. 

Work continued on model development during this month. Last month’s 
report discussed a simple code which was developed to check the validity of the model 
equations. In the simple code, the flow was assumed to experience constant accelera- 
tion as it moved through the nozzle. The code then was used to determine the variation 
of other properties across the nozzle; in particular, the area variation required for ob- 
taining the specified constant acceleration was obtained. 

The final code developed should, however, be capable of performing the 
reverse task. It should use a specified area distribution and determine the variations of 
velocity and other properties across the nozzle. This development was performed this 
month. To accomplish this, the model equations were cast in a modified form. Without 
this modification, solution of the equations would encounter numerical difficulties at the 
nozzle location where the gas speed became sonic. It should also be noted that the solu- 
tion could not be commenced at the nozzle entry. The gas velocity was zero at this loca- 
tion and, therefore, the area approached infinity. Thus, the solution was obtained for 
some small initial length of the nozzle assuming constant acceleration. Then the speci- 
fied area distribution was used to continue the integration to the nozzle exit. 

The code described above is currently being used to perform parametric 

studies. 


n., ir 

L.L. Narayanaswami r 



SUPERSONIC GAS-LIQUID CLEANING SYSTEM 


MONTHLY PROGRESS REPORT 

This progress report is submitted by Embry-Riddle Aeronautical Univer- 
sity (ERAU) to Bethune Cookman College (BCC), and describes the work done during 
the fifth month under TRDA Grant #41 1 towards the development of a Supersonic Gas- 
Liquid Cleaning System. 

Parametric studies were started this month using the model. The nozzle 
contour, see Figure 1, utilized in the cleaning system was used in these studies. 



Figure 1 : Nozzle Profile 


As the integration was performed, numerical difficulties were experienced near the 
location where the Mach number was 1 . The reason for this difficulty was that a variable 
Z (defined in the model) and its derivative dZ/dx* approached zero near the location 
where the Mach number was 1 . Owing to this problem, parametric studies were not 
completed this month. 


Currently, integration is being carried out using very small step sizes near 
the M = 1 location to overcome the problem. 


Last week (week of Feb 6, 1995 to Feb 10, 1995), Mr. Matt La Vigne of 
TRDA required that a write-up of the model work be sent to him. Accordingly, a report 
was prepared and submitted to TRDA. A copy of the report is attached. 


Also attached are summaries from some pertinent references located as 
a result of the literature search. 


LL. NarayanSswami 


r&b 





SUPERSONIC GAS-LIQUID CLEANING SYSTEM 
NOZZLE MODEL DESCRIPTION 


The model described in what follows describes a two-phase flow through a 
converging— diverging nozzle of a given geometry. The development closely follows that in 
Reference 1. As mentioned in the second monthly report, the model is based on the following 
assumptions: (1) the flow is steady and one— dimensional, (2) droplets are uniformly distrib- 
uted in the mixture, (3) droplets are uniform in diameter, and have uniform properties, (4) 
drag between the gas and liquid phases can be described in terms of a mean velocity differ- 
ence between them, (5) heat transfer between the phases can be related to the mean temper- 
ature difference between them and (6) there is no mass transfer between phases. 

Subject to these assumptions, the following equations can be used to describe 
the two phase flow: 


Mass ( Gas Phase ) : 

q u A = m 

(i) 

Mass ( Particle Phase ) : 

q p Up A = x m 

(2) 

Momentum ( Gas Phase) : 

+ 

8*1 •§• 
ii 

(3) 

Momentum ( Particle Phase) : 

■f 

EH#- 

II 

1 

(4) 

Energy ( Gas Phase) : Q 

dT dp 

u C p — u — = {u p - u) F p + Q p 

dx dx 

(5) 

Energy ( Particle Phase ) : 

Qp Up C s — Qp 

dx 

(6) 


In the above equations, Q and Qp represent the gas and particle densities, u and Up denote 
the gas and particle velocities, T and T p represent the gas and particle temperatures, and p 
is the pressure, m and xm denote the gas and particle mass flow rates, and Cp and C s are the 
gas and particle constant pressure specific heats. Utilizing assumptions (4) and (5), the drag 
force F p and the heat transfer Q p between the phases can be written as 


'HultbergJ.A. and Soo, SL..”Two-Phase Flow through a Nozzle,” Astronautica Acta, Vbl. 11. No.3, 1965. 



Drag : 


F p = f(u - u p ) 


(7) 


Heat Transfer : Q p = q (T — T p ) 


( 8 ) 


where f and q are constants related to the gas and particle phase properties. Equations (3) 
and (4) may now be combined to get the overall momentum equation as follows: 


du du p dp 

6 u — + q p Up + — — 0 

dx p P dx dx 


(9) 


Equations (5) and (6) could be used, in conjunction with Equations (3) and (4), to obtain the 
following overall energy equation: 


u 


Ur, 


(Cp T + — ) + X (Cs Tp + ) = constant 


( 10 ) 


In the above equation, the constant can be determined by using the nozzle boundary condi- 
tions. 


NONDIMENSIONAUZING THE MODEL EQUATIONS 


It is useful to nondimensionalize the equations governing the system. Accord- 
ingly, the various equations are nondimensionalized using reference variables po, To, Qo> u 0i 
L and (m/Qouo). In this list, po is the upstream tank pr essure, To is the upstream tank tempera- 
ture, 0o is the upstream tank density, uo is \/2CpTo, and L is the length of the nozzle. Thus, 
the equations are cast in terms of the following nondimensional variables: 


P 

Normalized pressure p* = “ 


(11 a) 



(lib)' 


T 

Normalized Gas Temperature T* = ZT 

+ fl 


Normalized Particle Temperature T p * 

II 

(11c) 

u 

Normalized Gas Velocity u * — 


(lid) 

Normalized Particle Velocity u p * = 

Up 

~“o 

(11c) 

A 

Normalized Area A * - 


(11/) 

X 

Normalized Length x * = ~j~ 


(11*) 


The normalized form of the equations are given below. 


Gas Mass : p* u* A* = T* 


* ( 12 ) 



du p * 



Particle Momentum : u p * 

— = 

C ( u * ~ 

Up*) 


dx* 



2 C p 

du* 

2 C p 

du p * 

Overall Momentum : u * 

— + 

XU 

* 

R T* 

dx* 

R T* 

dx * 


dT p * 




Particle Energy : u p * = D ( T* — T p *) 

dx* 


1 dp* 
P* dx* 


(13) 


= 0 (14) 


(15) 


Overall Energy : T * 


C s 

+ u* 2 + x — T p * + x u p 

C p 


C S 

*2 = 1 + x — 


(16) 



In the above set of equations, C and D are nondimensional equivalents of f and 
a respectively. Further, it should be noted that the ’constant’ on the right hand side of the 
dimensional form of the overall energy equation (Equauon (10)) ha s been obtained by ob- 
serving that, at the entrance to the nozzle from the tank (x* = 0), u and u p are zero, and 
T* an d T p * are 1. This results in the ’constant ’ being evaluated to be (1 + x C s /C p ) as shown 

in Equation (16). 

The five equations (12) through (16) describe the variation with x* of six vari- 
ables o* T* T * u* u n * and A*. If the variation of one of these variables with x* is known, 
then the five equations can be used to obtain the variation of the other five with respect to 
x*. In particular, the flow properties at the exit plane of the nozzle can be obtained. 


INITIAL calculations using model equations 

In order to verify the accuracy of the model equations, preliminary calculations 
were performed by assuming the gas velocity to be known through the nozzle. A simple varia- 
tion namely that u* varies linearly with x* was assumed (u* = ox^. Then, Equauons (12) 
through (16) were used to determine the variation of the other variables with x . The ap- 
propriate^ boundary conditions at the nozzle entrance (x* = 0) were p* = 1 , and, as menuoned 
in connection with Equation (1 6) above, u p * = 0, T* = 1 and T p - 1. 

In performing these calculations, the foUowing values were ^sumed for the 

several constants appearing in the model: D * 10., C = 100., C p - 1004.5 J/kg K, C s 9 
T/kg-K R = 287 J/kg-K and x = 0.3. In addition, the constant relating u* to x (the quantity 
a in u*’= ax*) was mken to be 0.9. The results of the calculations are shown in Figure 1. 


SOLUTION FOR ARBITRARY AREA RATIOS 


While the solution described in the previous section is useful as a reference operating 
characteristic, the variation of the gas velocity through the nozzle is notreadily available m pra« lca 
situations. More usually, A* would be known as a function of x*. Huts, it t is 
a solution methodology that would utilize the model equations (Equations (12) though . (1M I » de 
tetmine the variation of u*. u p *. T*, T p * and p* with x* for any aibttraty nozzle profile (A (x )). 
The methodology is described in this section. 

A close inspection of the equation set shows that singularities are present in the ex- 
pressions for dp*/dx*, dT*/dx* and du*/dx*. The singularities occur at the locanon in the nozzle 
where the Mach number (based on the speed of sound in gas only) 

eouations have to be recast in a form suitable for solution. This requires the introduction of an auxil 
Sy variable Z defined to be 1/2(M 2 -1) 2 , where M is the Mach number based on the speed of sound 
in gas. The resulting equation set is presented below: 






(13) 


u B 


du p * 

dx* 


= C ’ (u* - u*) 


dip* 

Up * - D ( T* - T p *) 

dx * 


(15) 


dZ 

dx* 


— = 2M 2 ( 
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T* T* 2 
xC s dT p * 
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) [ ( - 2 xu p * ) — 

A* dx* R dx* 


C P dx* 


2C p u* 2 xC s dT p * 

] + (M 2 - 1) ■ [— — + 2 xu p * — ] 
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yRT* 2 c p dx* 


dx 


(17) 


Cs 

T* + u* 2 + x — T p * + x u p * 2 = 1 + x — (16) 

Cr 


C s 

'P 


Z = - (M 2 - l) 2 (9/? M 2 = 1 ? J2Z 

2 

yR 

u* 2 = T*M 2 — 

2Cn 


(18 a&b) 


*(19) 


The above set of six equations ((13), (15), (16), (17), (18 a) and (19)) describes the variation 
of 7 variables (u p *, T p *, Z, A*, T\ u* and M) with x*, the distance along the nozzle. If the nozzle 
profile information (A*(x*)) is known, the above set can be solved to obtain the variation of 
the other six variables with x* by making use of the boundary conditions at x* = 0: 


T* = T p * = 1; u* = M = Up* = 0; Z = 0.5 (20) 


It should be noted that in performing the integration, the negative sign in Equa- 
tion (18 b) applies to subsonic flow, and the positive sign to supersonic flow. Thus, as integra- 
tion is carried out from x* = 0, the negative sign is used until the location where M reaches 



unity. The applicable sign in Equation (18 b) is then positive as solution is extended to the 
supersonic regime. 

It should also be noted that the integration cannot be started at x* of zero, sincfe 
A* goes to infinity there. Consequendy, an alternate solution has to be assumed over some 
initial distance in the nozzle. In this work, the linear gas velocity variation (u*= ax*), discussed 
elsewhere, is assumed in the initial region (say, from x* = 0 to x* = 0.03). The solution is then 
continued using the specified nozzle profile. 

To get the correct solution, the value of a cannot be arbitrarily chosen since 
each value of a corresponds to a particular nozzle profile (recall the particular A* variation 
obtained with a = 0.9, see Figure 1) . Thus, the chosen value of a has to be checked for correct- 
ness. If found to be incompatible with the specified nozzle profile, it has to be suitably cor- 
rected. This checking and correcting is performed at the nozzle throat, as described below. 

It can be shown that the gas velocity u* at the throat is given by the following 
analytical expression: 


diip * ~ l 

u t * 2 = a m * 2 [ 1 +* — ] ( 21 ) 

du* 

where u t * and a m * are, respectively, the gas speed at rite throat and the speed of sound in 
the mixture. By comparing the values of u t * obtained from the numerical solution with that 
from Equation (21), the correctness of the choice of a is determined. If the velocities do not 
agree, then an a new is obtained as 


U( exact num 

= a 0 id + <5a; da = rj — 


where rj is a constant less than one, and x t * denotes the throat location. 

It should be pointed out that numerical integration is difficult near the region 
where the Mach number M based on the speed of sound in the gas becomes unity. Around 
this condition, both Z and dZ/dx* approach zero, and the solution should be carefully ad- 
vanced from one x* to the next. 

To check the solution technique for arbitrary area ratios, the nozzle profile 
computed with u* = ax* was used. This profile was used as input to see if the technique recov- 
ered the linear velocity variation, and the other variations exactly as depicted in Figure 1. The 
solution was obtained by taking a to be 0.9 (thus, no iteration was necessary here using Equa- 
tions (21) and (22) to obtain the correct a). Figure 2 presents the results. It can be seen that 
the method has successfully computed the values of the different variables. 








SUMMARY OF REFERENCES 


’A NUMERICAL SIMULATION OF INJECTION OF DROPLETS IN A COMPRESSIBLE FLOW’, 
Daniel, E., Larini, M„ Loraud, J-C., Porterie, B. 

A numerical simulation of the injection of water droplets into a 2D, unsteady, inviscid 
compressible nozzle flow is computed. A two fluid (Eulerian— Eulerian) system was obtained by a weak 
conservative balance analysis of the gas phase and the dispersed phase. Each phase was considered a 
continuum. Drag, heat and mass transfers between the two phases are taken into account. 

Previous works by Chang (Ref 1. and 2.) considered two phase (gas— particle) flow in two 
and three dimensions. Golafshani et al (Ref 3.) had previously considered treating the dispersed phase 
as a discreet one and a Eulerian— Lagrangian model was solved through the use of a finite— volume 
method. The drawback to using such a method is the required interpolations between the Lagrangian 
grid (droplets) and the Eulerian grid (gas). The Eulerian— Eulerian method requires 

The governing equations were written in the weak conservative form. The energy 
equation for the dispersed phase reduces to the saturation temperature which indicates that the heat flux 
from the gas to the droplet is completely used in vaporization. This means that break-up or coalescence 
of droplets is neglected. An energy balance around one droplet is used to determine the mass transfer 
from the droplet. Yuen and Chen (Ref 4.) showed there is a diminuition of heat transfer due to mass 
exchange. A correction to the classical correlation of die Nusselt number was used to obtain the 
expression of the heat transfer. This corrected Nusselt number correlation is also applied to the drag 
transfer (Ref 5.) All the thermophysics variables are assumed constant except the latent heat of 
vaporization, the saturation temperature and the droplet enthalpy, all of which depend on the gas 
pressure. 

The MacCormack finite difference scheme was used to solve the partial differential 
equations. The boundary conditions are detailed for both phases. 

The results were not compared to experimental data because no experimental data was 
available for the injection of droplets into compressible nozzle flow. 

The results for steady flow were qualitatively compared to previous two-phase nozzle 
flow works. The two phase nozzle flow compared favorably with that of Chang. The droplet size was 
found to be the most important factor in affecting the droplet free zone in the nozzle. The void fraction 
was found to have no effect on this zone. 

Contours of the normalized number of droplets per unit volume are plotted for several 
times after droplet injection until the steady state is reached. In this paper the object of the droplet 
injection was to cool the gas and gas temperature contours are plotted for a sample nozzle at several 
times until the steady state was reached. 
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1 . I-Shin Chang. 'One and two-phase flows’, AIAA Journal, Vol 18, 1980 

2. I-Shin Chang. ’Three dimensional two-phase supersonic nozzle flows’, AIAA Journal, 

Vol 21, 1983 

3 . Mehdi Golafsham, Hai-Tien Loh, ’Computation of two-phase viscous flow in solid 
rocket motors using a flux split Eulerian— Lagrangian Technique*, aiaa/asme/asee 25th 
Joint Propulsion Conference, Monterey California 

4 m. C. Yuen, L. W. Chen, ’Heat transfer measurements of evaporating liquid droplets' , 
International Journal of Heat and Mass Transfer, Vol 27, 1977 
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transfer’. International Journal of Heat and Mass Transfer, Vol 27, 1983 



‘ONE AND TWO PHASE NOZZLE FLOWS’, Chang, I-Shing., AIAA journal, Vbl. 18. No 12, 
December 1980 


Modem propulsion systems usually have small throat radii, steep wall gradients at the 
entrance or submerged configuration nozzles in the interest of minimizing length and weight. A review 
of gas-particle n ozzl e flowfield investigations before 1962 was presented in Ref 2. More recent work 
includes the numerical iteration relaxation technique (Ref. 3) and an uncoupled flow model (Ref. 4). 
These studies help to explain some of the physical processes of gas-panicle transonic flows but all used 
the assumption of the fixed gas— phase streamline coordinates. This assumption is not true for steep 
gfpm«rif»s or small throat curvatures where the presence of particles will change the flow behavior. 
This paper analyses the time dependant, gas only and two-phase, inviscid, compressible flow in nozzles 
of arbitrary geometry. 

The equations of state are written in the weak conservative form for the unsteady, 
invoiced, transonic gas only and gas-particle flow inside nozzles with small throat radii of curvature, 
steep wall gradients and submerged configurations. The' friction, heat, drag, energy and momentum 
transfers are taken into account Henderson's correlation equation (Ref. 6) and Stokes law of drag 
coefficient for spheres in creeping motion are used in the momentum transfer. The particle Nusselt 
number is used in the heat-transfer (Ref 5.). 

Successive over-relaxation was utilized to generate the boundary fitted coordinates, BFC, 
(Ref. 7) through the use of the TOMCAT computer program. The FATCAT program (Ref 7.) was 
miiiwt to determine the factors needed to transform from the physical grid to the computational 
plane. The MacCormack finite difference scheme (Ref 8.) was applied to solve the transformed 
differential equations. 

For one phase flow, the initial condition was based on a one-dimensional, is en tropic flow 
analysis. An interpolation was used to obtain the gas only fkwresults and these were used as the initial 
»ctimaf«» for the gas-phase in the two-phase flow. The two-phase boundary conditions are discussed. 
The fourth order damping terms to the second order MacCormack method were retained to prevent 
nonlinear instability (Ref. 9). 

The solution method was applied to JPL nozzle, the Titan in solid rocket motor and the 
submerged nozzle configuration in the inertial upper stage solid rocket motor. This led to the following 
applicable conclusions: the MacCormack finite-difference method results in stable integration for 'both 
one and two-phase nozzle flow equations; BFC system improves the program capability to allow 
solution of flows inside complex geometry nozzles; small size particles stow down the gas phase 
expansion more than large size particles for a constant particle mass fraction; a high particle loading 
ratio can the gas phase to be subsonic at the throat; a constant fractional lag is not be justified for a 
two-phase transonic flow and the requires the momentum and energy exchange to be taken into account. 
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’TWO PHASE FLOW THROUGH A NOZZLE’, Hultberg, J. A., Soo, S. L.. Astronautica ACTA, Vol 
II, No 3, 1965 

Gas-particle two phase, one dimensional, supersonic flow in a nozzle was analyzed for 
spherical particles. Drag force, momentum and heat transfer between the phases was considered. The 
throat condition was obtained by setting the change in area with respect to distance measured along the 
longitudinal axis equal to zero. This led to the throat condition that the gas velocity is subsonic at the 
throat even when wall friction is neglected. The Runge-Kutta-Gill method (Ref. 1) was used. The 
boundary conditions were determined. The theoretical results were compared to experimental data and 
were found to be si gnifican tly different The one dimensional flow assumption was shown to be 
inadequate when compared to the experimental data. 
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’THREE DIMENSIONAL, TWO-PHASE SUPERSONIC NOZZLE FLOWS’, I-Shin Chang. AIAA 
Journal, Vol 21, No 5, May 1983 

The governing equations for a gas-particle, three dimensional, steady flow were written in 
the weak conservative form. The exchange of friction, heat, energy and momentum between the two 
phases was taken into account. Once again Stokes drag coefficient and Henderson’s correlations were 
used (Ref. 1 and 2). 

The equations were transformed into a cylindrical coordinate system (Ref. 3) before 
solution through die use of the MacCormack method of finite differences (Ref. 4.) A stability analysis 
calculation procedure (Ref. 3) was utilized. A tangent condition was applied for both the gas and the 
particle phase at the boundary surface. The Abbett scheme (Ref. 5) for determining gas phase boundary 
flow variables was modified to account for entropy changes at the boundary. Further boundary flow 
conditions are discussed. 
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’SUPERSONIC GAS-LIQUID CLEANING SYSTEM’, Caimi, Raoul E. B.. Thaxton, Eric A. 

Introduces the supersonic gas-liquid cleaning concept to replace the current hot, solvent 
flush method used to clean liquid oxygen systems at Kennedy Space Center. The performance which 
must exceeded in order to replace the current system are detailed. A handheld prototype system was 
build and successfully tested. 

One dimensional, irrotational, frictionless gas flow equations were used with a dispersed 
phase of droplets. Homogeneous, adiabatic, uniform velocity flow was assumed. The equation of state 
and the expansion law assuming an i sen tropic process was derived. An isentropic flow table was 
constructed and this was used to determine the Mach number versus area ratio. 

A H andhe ld prototype system using nitrogen and distilled water built and successfully 
tested. Additional applications were discussed. 


’A COMPARISON BETWEEN CENTERED AND UPWIND SCHEMES FOR TWO-PHASE 
COMPRESSIBLE FLOWS’, Daniel, E., SaureL R.. Larini, M., Loraud, J. G, 

AIAA/SAE/ASME/ASEE 29th Joint Propulsion Conference, June 1993, Monterey, California 

Two different methods were developed to analyze two dimensional, two-phase flow 
through an Eulerian approach. The different assumptions used which lead to conservative or 
nonconservative systems and to hyperbolic or degenerative hyperbolic systems are discussed.. One 
flow system was solved by a central scheme and a second system was solved by an upwind scheme. A 
Riemann solver for a dilute dispersed was included in the upwind scheme. The results of the two 
different schemes were found to be comparable. The differences in the boundary conditions were 
thought to have led to the differences between the results. The-Tipwind scheme was shown to accurately 
solve two-phase PDE in complex geometries. Hie central scheme could solve any two-phase flow that 
can be described 
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’EFFECTS OF COMPRESSIBILITY ON ATOMIZATION OF LIQUID JETS’, Zhou, Z. W„ Lin, S. P.. 
AIAA 30th Aerospace Sciences Meeting & Exhibit, Reno, Nevada, 1992 

The effects of compressibility of fluids on the atomization of a liquid jet emanating from a 
nozzle into the ambient gas was investigate through the use of a linear stability analysis. Increased 
compressibility of the liquid jet was shown to result in smaller droplets generated at a lower rate. It was 
alcn shown that small changes in compressibility result in a qualitative change of the mechanism of 
atomization. 



'NUMERICAL METHODOLOGIES FOR THE COMPRESSIBLE NAVIER-STOKES EQUATIONS 
FOR TWO-PHASE FLOWS’, Grasso, F., Magi, V., Modem Research Topics in Aerospace Propulsion, 

Pringer-Verlag 

In this large paper (33 pages) the physics of sprays and numerical methods for the solution 
of the compressible, two-phase Navier-Stokes equations are reviewed and a new numerical method is 
derived. The results from this method arc compared to experimental data for sprays in selected 
conditions. 
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SUPERSONIC GAS-LIQUID CLEANING SYSTEM 


PROJECT STATUS REPORT \p, a 

JANUARY 12, 1995 

This project status report is submitted by Precision Fabricating and Cleaning Co., Inc. 
(PFC) to Betbune-Cookman College (B-CC), and describes work performed between December 
8, 1994 and January 10, 1995 in support of this research project (Ref. TRDA Grant #411). 

Dr. David and Tom Jezowskd have begun testing of the nozzle system. This testing is 
taking place at Precision Fabricating and Cleaning in their Mobile Clean Room Facility. 

On December 14, 1994 we received the Breathing Air tube bank trailer (GT-13) here at 
PFC. The NASA prototype lab finished fabrication of the second nozzle on December 16 and we 
picked h up. The hot-wire anemometer became available to us on December 19 but it did not 
come with any type of probe. Dr. David chose a hot-film type of probe and B-CC procured it 
from the manufacturer. The probe arrived at B-CC around the end of December. 

All the equipment setup and connections were completed by January 3, 1995 and Dr 
David brought the hot-film probe to PFC on January 4, 1995. We began testing and were able to 
get consistent, repeatable velocity measurements. So fiir testing has been performed using gas 
onfy, at pressures ranging from 200 to 500 psig. Additional testing was performed on January 6th 
and 10th. Testing will continue on January 12th and 13th 


Eiley and Kathleen Hater came out on January 6th to observe our equipment 
setup. Unfortunately we had temporary equipment problems and could not demonstrate full scale 

testing while they were here, but we did fix the problem and resume testing the following 
morning. 

On January 10th we finished fabrication of the water injection assembly and connected it 
Into the existing nozzle system. Dr. David and Tom Jezowski began initial testing using the gas- 
liquid mixture. Further testing will continue with this setup. 


The next phase of testing will incorporate the use of either a supersonic phot tube or a 
Laser Doppler Velocimeter or maybe both. This will be determined by availability of the 
equipment And whether or not it is necessary. 

The data gathered so for, and in the fiiture, will be used by Dr.’s David and 
Narayanaswami to verify and “fine tune” the computer model developed by Dr. Narayanaswand 

On« the model is folly developed it will be used to optimize the nozzle design prior to actual 
fabrication of another nozzle. 

Maria Littlefield, Kathleen Harer, Kristen Riley, Eric Thaxton, and possibly others are 

sceduled to come out to PFC on January 12 to observe die equipment setup and testing b eing 
performed by Dr. David and Tom Jezowski. 


Precision Fabricating and Cleaning Co., Inc. 
Page 1 of 1 
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SUPERSONIC GAS-LIQUID CLEANING SYSTEM 


PROJECT STATUS REPORT Nn « 

FEBRUARY 13, 1995 

l Ste ^ S r ? >0rt * s ” bmittcd ^ ^««>n Fabricating and Cleaning Co., 

i ^ ii °®f Un r5?° kman Co ^ e 8 e ( B -CC), describes woric performed between 
Grant*# 41 1)"* ^ I ebruar >' 10 > 1995 in support of this research project (Ref TRDA 

Initial testing of the gas-liquid mixture was begun on January JJ 1995 

Observation of the flow characteristics was the first stage of this testing. During these 

observations the hot-film anemometer probe was inadvertently destroyed. The probe was 

not repairable. The data gathered with the probe using gas only was quite consistent and 

£ tc< * A,S< J» now that Or. David has received the calibration curve from Dan tec 
this data can be validated. 


> c 


Eflotts were made to locate a supersonic pitot tube but these efforts were 

w!SS? SfuI ^ i? VC ? 8tt ^ 1x1X0 prpcurin & l yP® ot P«>^ proved to be cost 
prohibitive. Dr. David and Tom Jezowski fabricated a pitot tube at P.F.C.’s facility >iamo 

small copper tiibing and a hypodermic needle. This phot tube was coupled to a pressure 
transducer and strip chart recorder supplied by NASA. 

This pitot tube configuration was used to gather data from the gas-liquid stream. 
Numerous pressure readings were taken over several days and Dr. David is in the process 
of analyzing this data. The nozzle system, pitot tube, and data gathering equipment 

remain m place in case further testing is required to verify any of die data already 
gathered. J 


Several black and white photographs of the shock waves exiting the nozzle were 
taken. The photographs turned out well. Even with water introduced into the system the 
shock waves were still visible. 7 

T!\ e next phase of testing will involve the use of a LDV located at NASA Langley. 
Dr. David, Maria Littlefield, and Eric Thaxton have already made one trip up to Langley 
to discuss this type of testing with the personnel familiar witli this equipment l believe 
arrangements are currently being made to transport the nozzle system up to Langley. 

The data gathered from the experimentation perfoimcd at P.F.C. will be used to 
fine tune the computer model developed by Dr. Naiayanaswami. This model will then be 
used to determine the optimal geometric design of the nozzle. Once this optimal design is 
determined, this nozzle should be fabricated and tested. 


Precision Fabricating and Cleaning Co., Inc. Jr / a (1* t 
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ATTACHMENT “B” 


Supersonic Gas-Liquid Cleaning System Research Project 
Final Project Budget Summary 



Research Project Entitled: 


SUPERSONIC GAS-LIQUID CLEANING SYSTEM 


FINAL ACTUAL PROJECT BUDGET 


Project Period: August 01 , 1 994 - January 31 , 1 995 


A. Direct Labor 

TRDA Program Director (15% time) 

TRDA Contract Manager (15% time) 

TRDA Secretarial Support (5% time) 

B. Travel (between TRDA, B-CC, E-RAU, PFC) 

C. Other Direct Costs 

D. Indirect Costs (TRDA @5%) 

E. Subcontracts 

1 Bethune-Cookman College 

a. Contract Mgmt. 

b. Faculty Stipends (3) 

c. Co-Op Students (3) 

d. Travel 

e. Subcontracts 

1 . Embry-Riddle Aero. Univ. 

(assistance in theoretical 
analysis of system) 

2. Precision Fabricating & 
Cleaning Co., Inc. 

(assistance in commercial 
configurations/testing) 

f. Indirect Costs (B-CC) 

F. Equipment 

TOTAL COSTS: 


* Extended to September 11,1 995 


NASA 

SHARE 

TRDA 

SHARE 

TOTAL 

COST 

$5,000 

$12,453 

$17,453 

$0 

$56 

$56 

$0 

$2,898 

$2,898 

$0 

$0 

$0 

$15,866 


$15,866 

$12,000 


$12,000 

$1,480 


$1,480 

$2,123 


$2,123 

$25,000 


$25,000 

$20,000 


$20,000 

$1,614 


$1,614 

$815 


$815 

$83,898 

$15,407 

$99,306 


$ 1 6,1 02 unspent NASA funds 

returned (de-obli gated) to NASA 








Research Project Entitled: 


SUPERSONIC GAS-LIQUID CLEANING SYSTEM 


PROPOSED BUDGET 


Project Period: August 01 , 1 994 - January 31 , 1 995 


A. Direct Labor 

TRDA Program Director (15% time) 

TRDA Contract Manager (15% time) 

TRDA Secretarial Support (5% time) 

B. Travel (between TRDA, B-CC, E-RAU, PFC) 

C. Consultants (as needed - technology transfer, 

engineering, project consultants) 

D. Indirect Costs (TRDA @5%) 

E. Subcontracts 

1 Bethune-Cookman College 

a. Contract Mgmt. 

b. Faculty Stipends (3) 

c. Co-Op Students (3) 

d. Travel 

e. Subcontracts 

1 . Embry-Riddle Aero. Univ. 

(assistance in theoretical 
analysis of system) 

2. Precision Fabricating & 
Cleaning Co., Inc. 

(assistance in commercial 
configurations/testing) 

f. Indirect Costs (B-CC) 


TOTAL COSTS: 


* Extended to September 1 1, 1995 



NASA 

SHARE 

TRDA 

SHARE 

TOTAL 

COST 


$2,129 

$2,871 

$0 

$3,530 

$574 

$595 

$5,659 

$3,445 

$595 


$0 

$500 

$500 


$0 

$450 

$450 


$0 

$5,000 

$5,000 

$12,500 

$18,000 

$15,000 

$2,000 

$95,000 

$0 

$95,000 

$25,000 
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SUPER-SONIC NOZZLE PROJECT 


1. BACKGROUND 

A proof-of-concept prototype of the supersonic gas-liquid cleaning system has been 
developed 1 and demonstrated at NASA, Kennedy Space Center. Presently, this system is 
being implemented as a cleanliness verification tool in replacement of CFC-1 1 3 rinsing. The 
system uses the high momentum of a supersonic air-liquid jet to remove contaminants from 
components while simultaneously emulsifying those contaminant into liquid jet (typically 
water). The liquid can be collected and sampled for contaminants to verify cleanliness. 
Extremely low volumetric flow rates of liquid are required with this system. NASA looked 
for a company that is interested in developing this system for commercial cleaning and 
cleanliness verification operations. The schematic diagram of the system is given in Fig.1 
and the mechanical dimensions of the prototype nozzle are shown in Fig.2. 

Under a research project of Technology Research and Development Authority (TRDA), 
NASA/KSC, Precision Fabrication &Cleaning company, Bethune-Cookman College (B-CC), 
and Embry-Riddle Aeronautical University (E-RAU) has jointly conducted study consisting 
of analysis and testing of the prototype cleaning system. Various tasks performed by these 
members are as follows: 

1. Modeling of fluid mechanics and nozzle flows, 

2. Parametric studies using the model, 

3. Optimization of design, 

4. Application to different matrixes, 

5. Commercial nozzle configurations. 

Embry-Riddle Aeronautical University (ERAU) did the modeling of the fluid mechanics and 
nozzle flows. Under the assumptions of steady state one dimensional flow, uniform droplet 
properties, drag between liquid and gas phases, heat and mass transfers between the 
phases, equations for mass, momentum and energy can be written for the two phases. In 
the extreme case of the liquid and gas phases having the same velocity and temperature 
(zero lag), the equations used by Caimi and Thaxton 2 are recovered. If the lag is small, but 
not zero, perturbation techniques provide the needed solution to the problem. In general, 
however, the lag will not be zero or small, and numerical methods are required to solve the 
equations in this work, this general condition of finite (not zero or small) lag is assumed. 

Bethune-Cookman College (B-CC) designed and conducted experiments in collaboration 
with PF&C and NASA/Langley Research Center (NASA/LaRQ through NASA/KSC to obtain 
the required data to confirm the theoretical model(s) and the validity of the assumptions. 
In order to do this, measurement of the nozzle flow parameters were done at different 
values of the independent parameters. 

This technical report presents details of these experimental investigations. 
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Fig.1 Schematic diagram of the Supersonic Gas-Liquid Cleaning System 





Fig.2 Mechanical dimensions of the mouth of the proto-type nozzle 
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2. The Supersonic Nozzle 

The nozzle is the main part of the supersonic gas-liquid cleaning system(SGLCS). The 
geometrical dimensions of the prototype nozzle under investigation are given in Fig-2. This 
converging-diverging type nozzle is made of steel and has throat diameter of 0.0812". 
Normally the SGLCS is operated and with upstream pressure (tank pressure) of nozzle 
around 300 psig. It can be shown that velocity of a compressible fluids like air and air- 
water mixture (with very small amount of water) will be in supersonic range at the exit 
plane of the nozzle. 

The flow is two-phased flow of air-water mixture. In Fig. 3, shadowgraph of air and air- 
water flows from the nozzle can be seen. The qualitative observation of the flow can be 
described as follows: 

1. At exit plane temperature of the flow seem to fall suddenly to a very low temperature 
from its initial temperature which can be taken as tank temperature. While moving away 
from the nozzle, it rises to be in equilibrium with ambient temperature. 

2. No water droplets or ice-droplet are found if the jet stream of the nozzle is not 
interrupted by any object. 

3. When an obstacle is placed in the stream, more condensation seems to occur at an axial 
distance greater than 10 cm as compared to the smaller ones. 

4. Best working distance seems to be 30 to 80 millimeters from the exit plane. 


3. Literature Survey 

The Literature Survey 3 revealed that most of theoretical and experimental studies reported 
in various journal were for subsonic fluid flows of the two-phased flow which are of interest 
to many industrial applications. However, the Anemometry and particle sizing are routinely 
done around in few high-speed wind tunnel experiments in aeronautical/aerospace 
research. These measurements are normally involved highly complex, expensive, and non 
portable custom-made equipment. This include Laser Doppler Anemometers, specially 
designed pitot probes, etc. It was observed that the measurement of the flow parameters 
such as, speed, direction, particle size, turbulence, vorticity, temperature etc in supersonic 
range put theoretical restrictions to the assumptions under which these parameters are 
measured in subsonic range by the conventional instruments. 

In this project, we looked for low-cost equipment and experimental methods which could 
be easily available and provide us the required data within the given time frame of the 
project. Following sections of this technical report describes the methods which were 
adopted for the purpose. 
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Figure 3 Shadowgraphs of supersonic flow of the nozzle at 300 psig. The pitot probe 
and linear scale can also be seen, (a) for breathing air, (b) for air-water mixture. 
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4. Measurement of Supersonic Flow Velocity 

The velocity of the nozzle stream in supersonic regime can be measured simply by a pitot 
probe, hot-wire/hot-film anemometer probe, Laser Doppler AnemometerdDA), or by Laser 
transit Anemometery. In this project we have used pitot probe and Laser Transit 
Anemometer to measure supersonic velocity U x at an axial distance x from the exit plane. 

4.1 By Pitot Probe 


A pitot probe is a tube with a blunt end facing into the airstream. A typical pitot probe has 
an inside to outside diameter ratio of 1/2 to 3/4, and a length of 15 to 20 diameters that is 
aligned with the airstream. Pitot probes are easy to construct and gives quite reliable 
measurements. It is an absolute measurement device which does not require calibration. 
Though the device is intrusive to the fluid flow , it is considered suitable because only one 
dimensional velocity data is required for this project Fig. 4(a) shows the dimensions of the 
pitot probe which was designed and built at PF&C. 
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Figure 4 (a) dimensions of the Pitot Probe, (b) schematic diagram showing various flow parameters 


An open-ended tube facing into the stream always measures the stagnation pressure it sees. 
Above M - 1 .0 the shock wave that forms the ahead of the tube means that it sees not the 
free stream stagnation pressure but the stagnation pressure behind a normal shock. The new 
value is called pitot pressure . In Fig. 3 a typical shock wave can be seen in the 
shadowgraph of the flow when pitot probe was aligned with the stream for the 
measurement. 



6 


4.1.1 Calculation of the Mean Flow Velocity using Pitot Probe 


Using theory of compressible fluids it can be shown 4 that 
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Where p,j x , p t1x , T x , and T, are defined in Fig. 4 (b). This dimensionless ratio can be 
found in the Standard table 5 for the range 1 ^ M u ^10. 

Using this table value of M lx can be determined for each value of stagnation pressure ratio 
measured at axial distance x.. These values are listed in column 10 of the Table 1. 

Similarly, if T x is defined the temperature at a distance x from the exit plane and T, is the 
fluid temperature inside the tank, it can shown that 



by substituting the value of M 1x obtained by equation(l) into equation(2), temperature T x of 
the supersonic flow can be determined at the measured axial distances x. These values are 
presented in column 12 of the Table 1. 

In compessible fluid medium the speed of sound at an axial distance x is given by, 


a. - Ji K, T. 

where V - 1.4 

R, = 287.0 J/kg. K ( gas constant ) 


(3) 


Using equation(3) speed of sound a, can be found for every x set experimentally. These 
values are shown in column 13 of the Table 1. 


Finally the supersonic velocity of the flow can be calculated by equation(4) given below, 


U. - a..M, (4) 

These values are given in column 14 of the Table 1. Fig. 5 depicts the variation of the 
velocity of nozzle flow both for breathing air and air-water mixture. 

















































































































































Figure 5 Variation of The Flow U x of the Supersonic Nozzle 
with respect to axial distance x, 
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4.2 By Laser Transit Anemometry 

Laser Transit Anemometry(LTA) is based on the direct measurement of the transit times of 
seed particles in the flow as they cross the foci of a pair of laser beams (Fig. 6). The transit 
time of an individual particle, %l, in conjunction with the known beam separation, s, 
provides a measurement of the velocity magnitude, v, , of particle given by 


The beams are rotated about an axis equidistant from and parallel to the two beams, 
enabling various beam orientations to be set with respect to the flow field. 

To detect particles as they traverse the foci of the instrument, scattered light from individual 
particles is detected by two photomultiplier tubes(Fig. 7) connected to filter-discriminator 
/ pulse-center detector circuits. As particles pass through the sample volume, pulses are 
produced on two separate channels, one for each beam. The method used to extract transit 
time information consists of computing the discrete cross correlation function in real time, 
r 12 ( i ), between the input pulse streams occurring on the two channels, P,(t) and P 2 (t). This 
discrete correlation function is given by 

OO 

r„(X) “ 21 P,(t,)P 2 (t, +%) <« 

t-o 

and is represented by as correlogram with 256 "bins" (or channels) with each successive bin 
representing a delay time X - i ^VT. 

In practice, turbulent fluctuations of the flow vector at the measurement point results in the 
LTA missing the particles passing out side the acceptance angle defined by the beam 
diameters and separation. Therefore, LTA data acquisitions performed by orienting the two 
beams at up to fifteen different angles with respect to the estimated mean flow direction and 
computing the correlogram using equation (2) at each angle. In this way an ensemble of 
correlograms in the ( - C,<) domain representing slices through the flow field is generated 
and can be analyzed to yield the desired flow parameters. In the experiments conducted 
at Langley for this project, only one angle was considered which gave the best correlogram. 
( fiC - 241 .000). 

4.2.1 The LTA system 

The LTA system which was used for the project at Langley is a Model 104, of Spectron 
Development Laboratories, Inc. It is a microcomputer based system (Fig. 8) in which LTA 
system control, data acquisition and data processing is performed by dedicated 
microprocessor based sub-systems. The optical system, xyuv drive, prototype nozzle, 
focusing Laser beams can be seen in photograph of Fig. 9 and two foci of Argon ion laser 
beams can be seen in the photograph of Fig. 10. 
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Figure 6 LTA Sample Volume Geometry 
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Figure 7 LTA Optical Head Schematic 





Figure 8 Microcomputer bused I TA system 
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Figure 9 General Experimental set up at NASA/Langley Research Center 



Figure 10 Nozzle flow and foci of two Argon ion Laser beams 
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4.2.2 Results of the LTA Measurements 

A typical correlogram(one of the 53 computer outputs) is shown in Fig. 1 1 . The mean 
velocity of the flow versus distance for supersonic gas-liquid nozzle is shown in Fig. 12. 
Variation of mean velocity with respect to distance was observed for tank pressures 220, 
320, and 420 psig. 


5. Measurements of Mass Rates of Water Droplets 

The flow of air-water mixture is obtained by the prototype control panel shown 
schematically in Fig. 1 and actually in Fig. 9. A very small amount of water is continuously 
mixed with compressed air before emerging out from the nozzle under differential pressure 
of 20 psig. The differential pressures are obtained by separate orifices shown in Fig. 1 . 

To order to assist the calculation of percentage content of water droplets in the supersonic 
flow of air-water mixture following observations were made: 

1. At the tank pressure of 420 psig, 5 liters of water took 90 minutes to pass through the 
nozzle with uniform rate. This gives flow rate of water droplets 0.926 grams per second. / 

2. At the tank pressure of 320 psig , 4 liters of water took 77 minutes to pass through the 
nozzle with uniform rate. This gives flow rate of water droplets 0.866 grams per second. 


6. Visualization of the Flow 

Visualization of the flow was done by shadow graph method. The experimental 





Obtained Via Laser Transit Anemometer Measurements 
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7. Minimum size of the Water Droplets 

The LTA experiments at Langley gave a quasi-quantitative conclusion about the average size 
of water droplets in the supersonic flow of the nozzle. Considering the 'Mean Velocity 
versus Axial Distance curve for 320 psig tank pressure shown in Fig. 12. When 5 cc of 0.8 
micron size polystyrene latex (PSL) particles is mixed with the water, the curve with lesser 
velocity was obtained(dipicted by round dots). This drop in the velocity of flow indicates 
that particles are heavier than the in flow obtained without PSL. This means water droplet 
size is much smaller that 0.8 micron. 


8. Concluding Remarks 

Mean velocity of the nozzle remains almost constant up to axial distance of 50 mm (2 
inches). It was observed that increasing the tank pressure increases the mean velocity of the 
flow at the exit of the nozzle. The optimum tank pressure, however, seems to be around 
300 psig. Best working distance seems to be between 30 to 80 millimeters. The best 
operating distance can be confirmed by the theoretical model of the flow. 
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INTRODUCTION 


re P ort is submitted by Embry-Riddle Aeronautical University (ERAU) to 
Bethune Cookman College (BCC), and describes the work performed under TRDA Grant 
#411 towards the development of a Supersonic Gas-Liquid Cleaning System. Commercial 
cleaning systems currently in use utilize the chemical action of a solvent to remove contami- 
nants from the surface being cleaned. A vast majority of these systems utilize CFC-113 as 
the solvent. Due to its detrimental effect on the earth’s ozone layer, there is a major effort 
worldwide to eliminate the use of CFC-113 as much as possible. The development of the 
supersonic gas-liquid cleaning system is an effort in this direction. 


The supersonic gas liquid cleaning system is shown in schematic in Figure 1. 
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Figure 1: Schematic of Supersonic Gas-Liquid Cleaning System 

It consists of a pressurized gas supply and a water supply. Gas leaves its tank through two 
pipe ines. In one line, the gas passes through pressure-regulating valves and flows into a 
converging-diverging nozzle. The gas passing through the second line flows into the water 
tank, and forces the water into the pipeline. This results in the flow of a mixture of gas and 
water through the nozzle. Due to the drop in pressure through the nozzle, water is atomized 
mto droplets as it flows towards the nozzle exit. Furthermore, the momentum interaction be- 
tween the gases and the water results in the latter leaving the nozzle at a high speed The 

impingement of this high-velocity jet of ’’heavy” water droplets has the potential for cleaning 
surfaces by ’’impact.” 

. Limited testing of the system illustrated in Figure 1 was carried out at NASA- 

KSC. In one of these tests, the system was used on an already precision-cleaned surface with 
a view to flushing out and capturing the residual contaminants on it. The flush was then ana- 
yzed to determine the concentration of the residual contaminants. Results were a measure 
of the cleaning capability of the original precision cleaning process used on the surface. In 
other words, the gas-liquid cleaning system was used as a tool for the ’’verification” of the 
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effectiveness of the precision-cleaning process. The tests indicated that the system was a vi- 
able alternative to CFC-based ’’verification” methods. 

The success of the gas-liquid cleaning system as a ’’verification” tool naturally 
led to its consideration as a primary cleaning device. The first step in this process was to deter- 
mine conditions that would be optimum for cleaning purposes. Accordingly, a two-phase 
flow model was applied to the nozzle. 2 In this analysis, the two phases were assumed to be 
thermal equilibrium with one another throughout the nozzle. Additional modeling work was 
done at University of Florida. 3 

, The current program of work was initiated to extend the limited experimental 

and modeling work mentioned above. The work was divided into two segments BCC con- 
ducted experimental studies with the system, and ERAU was involved in the development 
of a theoretical model to characterize the flow through the nozzle. 


One of the earliest investigations of two-phase nozzle flows was carried out 
by Hultbeig and co-workers. 4 The investigators developed a model, and tested it by compar- 
ing the results with those from an experimental program specifically designed for the inves- 
tigation. Other investigators had also been involved in the study of two-phase nozzle flows. 
The work of Reference 5 was concerned with the behavior of large water droplets in high 
speed flow. Both theoretical and experimental studies were performed as part of this effort. 
The investigation of Reference 6 was concerned with the development of a model describing 
two-phase flow through a nozzle. The model considered droplet breakup, and momentum 
and energy transport between phases. The model was subsequently used to perform para- 
metric studies. The model results were tested against experimental data as well. One other 
two-phase investigation that could be mentioned here was the one performed by Marble. 7 

This work was directed at the study and optimization of two phase flows in rocket motor 
nozzles. 


All the four investigations cited in the last paragraph modeled the nozzle flow 
as one dimensional. More recently, several multi-dimensional models have been developed. 
Shih-Chang has developed an axisymmetric model 8 and a three-dimensional model? de- 
scribing two-phase nozzle flows. Other investigators have also developed multi-dimensional 
models. 10 - 11 


While it would have been worthwhile to develop a multi-dimensional model 
accounting for droplet breakup, the short duration of the project precluded such an effort. 
Thus, the model developed for this work was one-dimensional, and did not account for drop- 
let breakup. The subsequent sections of the report describe model development, methodolo- 
gy developed for solving the model equations, and comparison of the results with experimen- 
tal observations. 
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NOZZLE MODEL DESCRIPTION 


The model described in what follows describes a two-phase flow through a 
converging-diverging nozzle of a given geometry. The development closely follows that in 
Reference 4. The model is based on the following assumptions: (1) the flow is steady and 
one-dimensional, (2) droplets are uniformly distributed in the mixture, (3) droplets are uni- 
form in diameter, and have uniform properties, (4) drag between the gas and liquid phases 
can be described in terms of a mean velocity difference between them, (5) heat transfer be- 
tween the phases can be related to the mean temperature difference between them and (6) 
there is no mass transfer between phases. 

Subject to these assumptions, the following equations can be used to describe 
the two phase flow: 


Mass ( Gas Phase ) : 

q u A = m 

(1) 

Mass C Particle Phase) : 

Qp Up A = x m 

(2) 

Momentum ( Gas Phase ) : 

du dp 

Q u ~r + — = F p 
dx dx 

(3) 

Momentum (. Particle Phase) 

du p 

: Qp u p ~ = ~ F P 

dx 

(4) 

Energy ( Gas Phase) : 

dr dp 

QuCp dc ~ -<*-“>*> + Qp 

(5) 

Energy ( Particle Phase) : 

dT p 

Q P Up C w — = - Q p 

(6) 


In the above equations, q and Q p represent the densities of the gas- and particle-clouds, u 
and Up denote the gas and particle velocities, T and T p represent the gas and particle tempera- 
tures, and p is the pressure, m and xm denote the gas and particle mass flow rates, and C p 
and C w are the gas and particle constant pressure specific heats. Utilizing assumptions (4) 
and (5), the drag force F p and the heat transfer Q p between the phases can be written as 
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Drag : 


F P = f(u - u p ) 


(7) 


Heat Transfer : Q p = q (T - T p ) 


( 8 ) 


where f and q are constants related to the gas and particle phase properties. Equations (3) 
and (4) may now be combined to get the overall momentum equation as follows: 


du du p 

eU & +QpUp ~^ 


dp 

+ — = 0 

dx 


(9) 


Equations (5) and (6) could be used, in conjunction with Equations 0) and (4), to obtain the 
following overall energy equation: 


u 

C Cp T + —) 
2 


4' ft (C w T p 4- 



= constant 


( 10 ) 


In the above equation, the constant can be determined by using the nozzle boundary condi- 
tions. 


NONDIMENSIONALLZING THE MODEL EQUATIONS 


It is useful to nondimensionalize the equations governing the system. Accord- 
ingly, the various equations are nondimensionalized using reference variables p 0 , T 0 , q 0 , u 0) 
L and (m/Q 0 u 0 ). In this list, p 0 is the upstream tank pressure, T 0 is the upstream tank tempera- 
ture, Qo is the upstream tank density, u 0 is V 2C p T 0 , and L is the length of the nozzle. Thus, 
the equations are cast in terms of the following nondimensional variables: 


P 

Normalized pressure p * = — 

P o 


(Ha) 
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Normalized Gas Temperature T* = 


Normalized Particle Temperature T p * = 


u 


Normalized Gas Velocity u * = — 


Ur 


Normalized Particle Velocity u„* = 

p «o 


Normalized Area A* — 


Normalized Length x* = 


(m/e Q u 0 ) 

x 


The normalized form of the equations are given below. 


Gas Mass : p* u* A* — T* 


du p * 

Particle Momentum : u p * = C (u* — u p *) 

dx* p 

2 C p du* 2 C p du p * 1 dp* 

Overall Momentum : u * + xu * h — 

RT* dx* r t* dx* P* dx* 


dip* 

Particle Energy : u p * = D ( T* — T„ *) 

dx* p 


Overall Energy : T * + u * 2 


+ x 



T p * + x Up* 2 


+ X 


( lib) 

(11c) 

(lid) 

(He) 

( 11 /) 

( 11 *) 

( 12 ) 

(13) 

0 (14) 

(15) 

(16) 
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In the above set of equations, C and D are nondimensional equivalents of f and 
q, respectively, and are given by 4 - 7 


C 


9 

2 Q w o 2 u 0 ' 


D = 


3kL 

C^j 2 q w u 0 


( Yla&b ) 


where p and k are, respectively, the coefficients of viscosity and thermal conductivity of the 
gas phase, a is the particle radius, and q w is the density of the particle material. Further, it 
should be noted that the ’constant’ on the right hand side of the dimensional form of the over- 
all energy equation (Equation (10)) has been obtained by observing that, at the entrance to 
the nozzle from the tank (x* = 0), u* and u p * are zero, and T* and T p * are 1. This results in 
the ’constant ’ being evaluated to be (1 + x C s /C p ) as shown in Equation (16). 

The five equations (12) through (16) describe the variation with x* of six vari- 
ables p* T*, T p *, u* u p * and A*. If the variation of one of these variables with x* is known, 
then the five equations can be used to obtain the variation of the other five with respect to 
x* In particular, the flow properties at the exit plane of the nozzle can be obtained. 

INITIAL CALCULATIONS USING MODEL EQUATIONS 


In order to verify the accuracy of the model equations, preliminary calculations 
were performed by assuming the gas velocity to be known through the nozzle. A simple varia- 
tion, namely that u* varies linearly with x* was assumed (u* = ax*). Then, Equations (12) 
through (16) were used to determine the variation of the other variables (T* T p * u p * and A*) 
with x*. The appropriate boundary conditions at the nozzle entrance (x* = 0) were p* = 1, 
and, as mentioned in connection with Equation (16) above, u p * = 0, 1* = 1 and T p * = 1. 

In performing these calculations, the following values were assumed for the 
several constants appearing in the model: D = 10., C = 100., C p = 1004.5 J/kg-K, C w = 494 
J/kg-K, R = 287 J/kg-K and x = 0.3. In addition, the constant relating u* to x* (the quantity 
a in u* = ax*) was taken to be 0.9. The results of the calculations are shown in Figure 2. 

SOLUTION FOR ARBITRARY AREA RATIOS 

While the solution described in the previous section is useful as a reference 
operating characteristic, the variation of the gas velocity through the nozzle is not readily 
available in practical situations. More usually, A* would be known as a function of x*. Thus, 
it is necessary to develop a solution methodology that would utilize the model equations 
(Equations (12) through (16)) to determine the variation of u*, u p *, T* T p * and p* with x* for 
any arbitrary nozzle profile (A*(x*)). The methodology is described in this section. 

As a first step in obtaining the solution for arbitrary area ratios, the equations 
are written in a suitable form. The resulting equation set is presented below: 
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(13) 


du p * 

Up * = C (u* — Up*) 

dx* 


dT p * 
dx * 


= D (T* - T p *) 


(15) 


du * 
dx* 


T* dA* 2CpXu * 
A* dx* + ^ r 


2 xu p * ) 


du p * 

dx* 


T* 2 C p 

— - u* + 2u* 

U* R 


xC w dT p * 
C p dx* 


(18) 


T* + u * 2 


+ 



T p * + x Up * 2 


1 + x 



(16) 


The above set of four equations ((13), (15), (16) and (18)) describes the variation of 5 variables 
( U P*» Tp* , A*, T* and u*) with x*, the distance along the nozzle. If the nozzle profile informa- 
tion (A*(x*)) is known, the above set, in principle, can be solved to obtain the variation of 
the other four variables with x* by making use of the boundary conditions at x* = 0: 


T* = T p * = 1; u* = M = u p * = 0 


(19) 


Unfortunately, the integration cannot be started at x* of zero, since A* goes to 
infinity there. Consequently, an alternate solution has to be assumed over some initial dis- 
tance in the nozzle. In this work, the linear gas velocity variation (u* = ox*), discussed else- 
where, is assumed over some initial nozzle length (say, from x* = 0 to x* = 0.03). This process 
would provide a set of boundary conditions at x = 0.03. The solution is then continued using 
the specified nozzle profile. 

The value of a cannot, however, be arbitrarily chosen since the earlier discus- 
sion with linear velocity profiles shows that a particular value of a corresponds to a particular 
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nozzle profile (recall the particular A* variation obtained with a = 0.9, Figure 2). Thus, the 
value of a has to be checked for correctness. If found to be incompatible with the specified 
nozzle profile, it has to be suitably corrected. This checking and correcting is performed at 
the nozzle throat, as described below. 

It can be shown that the gas velocity u t * at the nozzle throat is given by the 
following analytical expression: 4 " J 


du p * -l 


u t * 2 = a m * 2 [ 1 + x ] 

du* 


( 20 ) 


where a m * is the speed of sound in the mixture. This expression is obtained by noting that 
dA*/dx* is zero at the throat location of the nozzle. Thus, there are two (one from numerical 
solution and the other from Equation (20)) independent estimates of u* at the nozzle throat 
location. The two values will agree if a is correcdy chosen. If they do not, a new value of 
a, a n ew, is obtained as 


a new &old = 1 1 


*Eqn.2Q *num 


X,’ 


( 21 ) 


where r\ is a constant, and x t * denotes the throat location. 


A computer code in FORTRAN has been developed to carry out the solution 
method oudined above. The code utilizes an Ordinary Differential Equation Solver to carry 
out the necessary integration. 


The solution technique for arbitrary area ratios was checked by applying it to 
the nozzle profile obtained by using the linear velocity variation. Accordingly, several linear 
velocity profiles were assumed (a = 0.6, a = 0.7, a = 0.8, a = 0.9). For each of these cases, 
the nozzle area variation was obtained. These area variations were then used as inputs to the 
arbitrary area ratio method to determine if it predicted back the linear velocity variations. 


In each of the cases tried, the method successfully predicted the correct linear 
velocity variations associated with the specified nozzle profiles. Table 1 below compares the 
values of a obtained from this method with their correct values. 


TABLE 1 : Comparison of and a num 


ttact 

a num 

0.9 

0.901 

0.8 

0.798 

0.7 

0.695 

0.6 

0.599 
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It should be noted that, with linear velocity variation in the nozzle, the check- 
condition (Equation 20) can be shown to reduce to 


( i *2 

i P 

1 + X- 

a 


with a m * 2 = ^-(l + a 2 x* 2 ^-) , 

^ eq 


P = UJC 2 + 4Ca - C ), and 


C«t = C P [D( 1 + ^) + 2fi](l + *^~\d + ^) _1 


( 22 ) 


(23) 

(24) 

(25) 


Since the actual nozzle velocity variation is known to be linear in the test cases, Equadon (22) 

has been used to determine the correct value of u t * required for updating the value of a (Equa- 
don (21)). 


It is appropriate, at this juncture, to mendon another solution method tried dur- 
ing the course of this investigation. In this method, u* was replaced by an auxiliary variable 
Z, defined as 1/2(M^— 1)^> where M is the Mach number based on the speed of sound in gas. 
The variation of Z with x* is described by the following set of equations: 


dZ 

dx* 


„ 1 u* 2 T* dA* 2 Cpxu 

= 2M 2 ( — + ) [ ( _ 

T * j*2 A* dx* 


- 2 xu p * ) 


dx 


xC w dT p * 2 C p u * 2 xC w dT 0 * 

+ — ] + (M 2 - 1 ) — 


du p * 


dx- 


yRT * 2 


w ui p —r 

[ 7 r +2 xu p * ] (26) 

c - dx* -*-* 


dx 
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OR M 2 = 1 =F JlZ 


(27 a&b) 


Z 



2 


yR 

u* 2 = T * m2 ^T (28) 

In Equation (27 b) above, the negative sign applies when flow velocities are subsonic, and 
the positive sign applies for supersonic flow velocities. 

Equations (13), (15), (16), (26), (27 a) and (28) are 6 equations describing the 
variation of 7 variables (u p * T p * Z, A* T*. u* and M) with x* If the nozzle profile (A*(x*)) 
is specified, then these equations can be solved, beginning with the following boundary conditions: 


T* * T p * = 1; u* = M = u p * = 0; Z = 0.5 (29) 

The integration performed with this set of equations encountered severe numerical difficul- 
ties near the location where the Mach number based on the speed of sound in gas approached 
unity. Several remedies, such as (1) the use of analytical jacobian in integration and (2) the 
use of L’Hospital’s rule, were tried but the problem persisted. This approach was, consequent- 
ly, abandoned, and the approach, described earlier, directly using the flow velocity u* was 
utilized to obtain solutions. 
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NUMERICAL SIMULATION OF NASA-KSC SUPERSONIC NOZZLE 

The solution procedure for arbitrary area variations, developed and verified 
above, was utilized to numerically simulate the flow behavior in the NASA-KSC supersonic 
nozzle. These efforts and the results therefrom are presented in this section. 



0.397” dia. 0.0812” dia. 0.1895” dia 



e 


Figure 3: Schematic of NASA-KSC Nozzle 


Figure 3 shows the NASA-KSC supersonic nozzle in schematic. As illustrated 
in, and described in connection with, Figure 1, the nozzle is supplied with a mixture of air 
and water. The diameter of the nozzle passage linearly decreases up to the throat, and then 
linearly increases. By maintaining a sufficiently high pressure upstream of the nozzle, the gas 
flow can be accelerated to supersonic speeds in the nozzle. It is expected that the gas would, 
by momentum interaction, accelerate the water droplets. Thus, under the right conditions! 
the exit jet from the nozzle consists of gas flow with entrained high-velocity water droplets, 
and may be used as a surface cleaning device that works through ’’impact.” 

The area profile of the above nozzle was used as input into the developed com- 
puter code in an effort to simulate the flow through it. Again, a guess value for a was used 
to start the solution (recall that the guess is to be checked, and, perhaps corrected, at the 
throat). However, no solution could be obtained. Initially the difficulty was thought to be due 
to the guess value for a being too far removed from its correct value. Consequently, a wide 
range of values for a was tried. Again, no solution could be found. After some reflection, the 
reason for this behavior was determined to be the incompatibility between the area and the 
velocity information specified at the upstream end of the nozzle. 

Consider a conveiging-diverging nozzle attached to an upstream tank. The 
passage area of the nozzle is initially very large (theoretically infinity at the junction between 
tank and nozzle). In the NASA-KSC nozzle, the diameter at nozzle entry is 0.397”. This corre- 
sponds to an area rado (area/throat area) of 23.9. The Mach number corresponding to this 
area rado is about 0.02, which corresponds to a velocity of about 7 m/s. While this is small, 
it is not small in the ’’numerical sense.” In carrying out the soludon for the NASA-KSC nozzle, 
the velocity at the nozzle entry is specified to be zero, and the area rado there is 23.9. These 
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two items of information ate fundamentally incompatible with one another. It is this funda- 
mental incompatibility that is responsible for the failure of the numerical method in this case. 
If the area variation were known from a ’’truly” low-velocity location, the numerical method 
would have been successful. Support for this statement is provided by the success of the 
method in obtaining the correct velocity profiles in the test cases. In each of these test cases, 
the area variations were known from a ’’truly” low-velocity location unlike in the NASA-KSC 
nozzle. 


Due to this difficulty, the method and the associated code, described above, 
are not of utility in numerically simulating the flow through the NASA-KSC nozzle. Conse- 
quendy, an alternate method has to be developed for this purpose. The alternate method, 
with the rationale behind it, is described in what follows. 


Since the flow is modeled as one-dimensional, the characteristics of the flow 
at any nozzle location depend on the area there, and not on the actual nozzle profile used 
to reach that location. Consider two nozzles both with the same exit-to-throat area ratio, but 
which use different profiles to arrive from the throat to the exit. Within the framework of one- 
dimensional flow, the exit conditions from both nozzles would be the same, in spite of the 
different nozzle profiles utilized in them. This observation suggests a straightforward method 
for numerically simulating flow through a given nozzle. 

In the numerical simulation, the velocity variation through the nozzle is as- 
sumed to be linear, that is u* at a distance x* from the nozzle entrance is assumed to be ax*. 
This assumed u* variation can be used in Equations (12) through (16) to obtain closed form 
solutions for the other significant variables (u p * T p * T*, p* and A*) as follows. 


= fix* 

(30) 

= 1 _D[ <> 2 + *0 2 I 2 

P 1 *c„ J* 

(31) 

+ — ) + 2 J 


= 1 -a 2 ** 2 ^. 

^ eq 

(32) 

= (l 

'-eq / 

(33) 

1 

(34) 

* /l 7 *2 CpOrd + *i) 

ax* (1 - a 2 ** 2 ^) 
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where P and C eq are given by Equations 24 and 25, respectively. Thus, it is possible to deter- 
mine the variation of u*, u p * T p * T* p* and A* with location x* along the nozzle. These equa- 
tions form the basis for a code that has been developed in FORTRAN for the numerical simula- 
tion of the supersonic nozzle. The code is named ssl.f. 

In performing the calculations, the actual value chosen for a has no effect on 
the flow properties at a given area ratio location. In fact, the program user may convince him- 
self of this behavior by running the code for two different values of a. Comparison of the 
results would show that the flow properties at locations with the same area ratios would be 
the same for the two cases. However, the x* locations for attainment of a given area ratio 
would not be the same in the two cases. A larger a would simulate a nozzle with more rapid 
acceleration than a smaller a, and thus would have a more rapidly varying area ratio with 
distance. Consequently, a given area ratio would be attained at a smaller x* in a larger a nozzle 
than in a smaller a nozzle. 

A study of the equations and the program will reveal that the results are inde- 
pendent of the pressure upstream of the nozzle. Specifically, the method assumes that the 
nozzle expands the flow to supersonic speeds, regardless of the pressure upstream of the 
nozzle. This is simply a result of the physical phenomenon that, once the upstream pressure 
is large enough to establish supersonic flow through the entire nozzle, the flow becomes in- 
dependent of the upstream pressure, and depends only on the nozzle area ratio. Thus, there 
exists a critical pressure above which the nozzle solution becomes independent of the up- 
stream pressure, and supersonic flow would always be maintained. A tacit assumption made 
in the present development is that the code will be exercised only above this critical pressure. 
This is a reasonable assumption to make since below the critical pressure, shocks enter the 
nozzle, and the resulting flow is not expected to be effective in cleaning due to the drop in 
flow velocities across shocks. 

The above-mentioned critical pressure is a function of the area ratio of the 
nozzle used. The critical pressure determination for a given area ratio nozzle is a straightfor- 
ward problem in gas dynamics. A FORTRAN computer code named pcrit.f has been devel- 
oped for this purpose, and accompanies the report. The code requires as input the value of 
the nozzle area ratio, and provides as output the minimum suggested pressure upstream of 
the nozzle. 


In the interest of completeness, a related problem should also be addressed 
at this point. Just as a minimum pressure exists for a given area-ratio nozzle, there exists a 
threshold area ratio that must not be exceeded for a given upstream pressure. If this area ratio 
is exceeded, shocks again would enter the nozzle. For operating pressures currently in use 
with the cleaning system (around 300 psig), this area ratio is very large, and need not be cause 
for concern. However, another computer code named acrit.f has been developed for deter- 
mining this maximum area ratio. The program requires as input the upstream pressure, and 
provides the threshold area ratio on output. 
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PARAMETRIC STUDIES WITH THE NASA-KSC NOZZLE 


The developed computer code was used to perform very limited parametric 
studies. It would have been desirable to carry out extensive investigations with the developed 
model. However, partly due to the short duration of the project, and partly due to the time 
consumed in developing the model, parametric studies were restricted to the two cases for 
which experimental comparisons were available. 

Accordingly, the NASA-KSC nozzle was chosen as the candidate for the para- 
metric studies. It was 0.5” long, and had throat- and exit-section diameters of 0.0812” and 
0.1895”, respectively. Air and water were fed into the nozzle from upstream tanks. The various 
controlling variables were chosen to be the same as those used in the experiments. The ex- 
periments were performed by Dr. Sunil David of BCC, and the reader is referred to Dr. David’s 
report for details of the experimental setup and measurements. As mentioned earlier, two 
experimental conditions were chosen for comparison with the model. In the first, the pressure 
upstream of the nozzle was 320 psig, and the air supply tank was at a temperature of 11°C. 
The water flow rate in this experiment was 0.0519 liters/min (4 liters in 77 minutes). In the 
second, the pressure upstream of the nozzle was 420 psig. The air tank temperature was again 
11°C. The water flow rate into the nozzle was 0.056 liters/min (5 liters in 90 minutes). The 
water flow rates in both cases provided very small particle-to gas mass loading (x of 0.0312 
in the first case and 0.0474 in the second). These values were input into the developed code 
to determine the nozzle flow characteristics. Of particular interest were the gas velocities at 
the nozzle exit. In constructing the solution, a was assumed to be 0.87 in both these cases. 

In addition, air was assumed to have a coefficient of viscosity p, of 1.7894 (10 -5 ) 
kg/m/s and a coefficient of thermal conductivity of 2.4335 (10 -6 ) J/m/s/K. The specific heats 
of air and water were taken to be 1004.5 J/kg-K and 4l84j/kg-K, respectively. The density 
of water was assumed to be 1000 kg/m 3 . Water droplets were assumed to be spherical, and 
each water droplet was assumed to have a diameter of 2.5 microns. Based on these assumed 
properties, the nondimensional coefficients C and D, defined in equations (17 a&b) were de- 
termined to be 0.87 and 1.88(10-5), respectively. 

Results from the code are presented in Figures 4 and 5. Figure 4 corresponds 
to the first case, and Figure 5 to the second. In each of these Figures is shown the value of 
x* where the area ratio is equal to the exit-to-throat area ratio of the NASA-KSC nozzle. Thus, 
the computed values of velocities at these locations should compare with the exit velocities 
from the NASA-KSC nozzle. In both cases, the computed air velocity is 670 m/s. The indepen- 
dence of the nozzle exit conditions to upstream pressure has been explained elsewhere in 
the report. The measured values are about the same, at 630 m/s, see Dr. David’s report, for 
the two cases. Considering the simplicity of the model, the excellent agreement between the 
model and experimental values is noteworthy. 

It can also be seen from the Figures that the air temperature rapidly drops 
through the nozzle in both cases. In contrast, the particle-phase temperature remains 
constant at its upstream value all through the nozzle. This result seems to be supported by 
experimental observations: Dr David has observed that the droplets issuing from the nozzle 
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are in liquid form, suggesting that they are at temperatures above 0°C. This result suggests 
that the initial model of Caimi and Thaxton 2 may not properly describe the nozzle flow, as 
explained below. 

In their development, Caimi and Thaxton 2 assumed thermal equilibrium (due 
to rapid heat transfer) between the two phases. Results of the current investigation would 
suggest that the assumption of no heat transfer between the phases would be much closer 
to reality than that of rapid heat transfer. 
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CONCLUSIONS 


This project was concerned with the investigation of a supersonic gas-liquid 
cleaning system. The work consisted of an experimental study of the flow of air-water mix- 
ture through a nozzle. This experimental work was supplemented by the development of 
a theoretical model capable of simulating the nozzle flow. Bethune Cookman College con- 
ducted the experimental part of the study, and ERAU conducted the theoretical part. 

The work has resulted in the successful development of a theoretical nozzle 
flow model. The associated computer code named ssl.f is provided with this report. It is easy 
to use, and there is a good possibility that it can be executed on a personal computer. In 
addition, two other codes pcrit.f and acrit.f have been included. The need for, and the use 
of, these codes have been explained in the report. Both of these codes utilize a nonlinear 
algebraic equation solver from the IMSL library. If the host computer does not have IMSL, 
then it would not be possible to recompile these two source codes to generate the required 
executables. To overcome this problem, in addition to the source program, the executables 
are also being supplied with this report. 

There are two diskettes attached. Each of these diskettes contains the same 
files; namely, ssl.f, pcrit.f, acrit.f, ssle, pcrite and acrite. ssl.f, pcrit.f and acrit.f are source codes, 
and ssle, pcrite and acrite are associated executables. One diskette contains the files in UNIX 
format, and the other in DOS format. The DOS— format diskette contains the six files, stored 
in succession. In the UNIX— format, the files have been gathered into a master tape archive 
file called ’exe'. This master archive contains two files, one named ’exe’ and the other named 
’source’. File ’exe’ contains the executable files, and file ’source’ contains the source codes. 
The following procedures may be used to retrieve and use these files in a UNIX-driven 
workstation. 

RETRTVING FILES FROM DISKETTE ONTO COMPUTER 

(1) Create a directory called nasa: mkdir nasa 

(2) Make ’pulse’ the working directory: cd nasa 

(3) Extract file ’exe’ and ’source’ from diskette: tar xvf /dev/rfdO exe source. There should 
now be two files named 'exe' and ’source’ in the working directory. 

(4) Extract files ssl.f, pcrit.f and ’acrit.f from ’source’: arx source. The three source files 
should now be in the directory. 

(5) Extract files ’ssle’, pcirte’ and ’acrite’ from ’exe’: ar x exe. The three executables should 
now be in the directory. 

COMPILING AND LINKING TO CREATE EXECUTABLES 

(1) Compile ’ssl.f and store object code in ’ssl.o’: f77 -c ssl.f 

(2) Compile ’pcrit.f and store object code in ’pcrit.o’: f77 -c pcrit.f 

(3) Compile ’acrit.f and store object code in ’acrit.o’: f77 -c acrit.f 

(4) Link to create executable ’ssle’: f77 -o ssle ssLo 

(5) Link to create executable ’pcrite’: f77 -o pcrite pcrit.o ’path’/IibimsLa. As mentioned 
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elsewhere, ’pcrit.f utilizes the IMSL library for solving a nonlinear equation that arises in the 
problem. Thus, if IMSL is not availabe in the host computer, this command would not work. 
If IMSL does reside in the host computer, the path leading to the directory holding the IMSL 
library is to be used in place of ’path’ in the command. 

(6) Link to create executable ’acrite’: f77 -o acrite acrit.o ’path’/liblmsl.a. Again, since this 
program uses IMSL, comments made in reference to creating ’pcdte’ apply here as well. 

EXECUTING CODES 

(1) To run ’ssle’: ssle. Execution creates two output files ’outl.dat’ and ’out2.dat’. 

(2) To run ’pcrite’: pcrite. Results are displayed on screen. 

(3) To run ’acrite’: acrite. Results are displayed on screen. 

Continuing with the concluding remarks, the nozzle model was exercised at the same condi- 
tions for which experimental data were obtained in order to determine its validity. It was 
found that the model correctly predicted the experimental behavior at the two conditions 
checked. The results also pointed out that the assumption, as made in Reference 2, of thermal 
equilibrium between the two phases may not be valid. 

It would have been desirable to conduct extensive parametric studies with the 
model. Time constraints did not permit such studies. Time constraints limited the scope of 
the experimental investigation as well. To more fully characterize the nozzle and its capability 
as a cleaning tool, ERAU may continue to work on this problem in the coming months. The 
work, if pursued, will likely employ the services of a Masters-level student providing good 
educational experience to the student on a project of considerable current interest. 
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SUPERSONIC GAS-LIQUID CLEANING SYSTEM 


AL PROJECT REPORT 


MARCH 22, 1995 


This project status report is submitted by Precision Fabricating and Cleaning Co., 
Inc. (PFC) to Bethune-Cookman College (B-CC), and summarizes work performed 
between September 12, 1994 and March 1 1, 1995 in support of this research project (Ref. 
TRDA Grant #411). 

This Supersonic Gas-Liquid Cleaning System research project came about in 
support of NASA’s technology transfer program. This approach to miltin g NA S A 
technology with local manufacturing reflects NASA’s involvement in programs and 
projects that are mutually beneficial to the parties involved. 

The development of this Supersonic Gas-Liquid Cleaning System came about in 
response to die mandatory elimination of CFC’s from cleaning operations. The CFC’s 
used in these operations were deemed to be environmentally hazardous and therefore 
alternative methods of cleaning and subsequent validation of cleaned items were required. 
NASA engineers developed this nozzle system as one alternative method of cleaning and 
validation. In fact, it’s use has already been approved in NASA KSC’s cleaning 
specification KSC-C-123G, “Surface Cleanliness of Fluid Systems, Specification for”, 
dated August 4, 1994. 

Even though it’s use was already approved in KSC-C-123G, the supersonic nozzle 
system still required a characterization study to verify the engineer’s claims of 
“supersonic” properties. Theoretically these supersonic properties exist but thus far no 
anpirical date had been gathered to substantiate these claims. Therefore one objective of 
the characterization study was to provide empirical verification of the theoretical design. 
In addition, only one nozzle configuration had ever been fabricated. Any information 

about two-phase flow through a nozzle that could improve the design of future nozzles 
was certainly desirable. 
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A second objective of the study was to determine an optimal geometric design of 
the nozzle for use in cleaning operations currently taking place at Kennedy Space Center. 
Whatever the application of die nozzle, its optimal design should call for the water 
droplets to possess the maximum velocity practically possible. That way the kinetic 
energy used for cleaning will come from die velocity of the droplets and will not be so 
dependent upon the mass of the droplets. The current nozzle configuration uses 
approximately one quart of water per one-half hour of operation at 300 psig. The optimal 
operating range for cleaning surfaces is currently approximately two inches (axially) from 
the exit plane of the nozzle. 

A third objective of this phase of die project was the development of a computer 
model in which system parameters could be varied to determine different optimal designs 
depending on application criteria. Inherent to this characterization study was die 
requirement for a literature search to locate any information that described current surface 
cleaning methods, information that would be of assistance in die optimization of nozzle 
designs, as well as any computer modeling that might have already been completed 
regarding two-phase flow through a nozzle. 

The accomplishment of the above mentioned three project objectives required the 
participation of several parties. Bethune-Cookman College (B-CC) and Embiy-Riddle 
Aeronautical University (ERAU) provided support with analytical verification, 
mathematical modeling, as well as experiment design. Precision Fabricating and 
Cleaning Company, Inc. (PFC) provided testing facilities, manpower, and system 
components to facilitate testing of a second nozzle system. NASA Kennedy Space Center 
(KSC) fabricated and supplied the actual second nozzle and water injector assembly 
required for this second nozzle system. NASA’s Langley Research Center (LARC) 
supplied facilities and manpower to support final stages of nozalc testing. 

First the nozzle testing procedures and equipment requirements were discussed 
and decided upon by doctorate professors from Bethune-Cookman College and Embiy- 
Riddle Aeronautical University. Analytical verification of the existing prototype nozzle 
system was earned out by ERAU. This involved verifying the equations, calculations, 
and investigating any assumptions made by the NASA design engineers. Also 
investigated was the feasibility of applying basic gas flow equations to two-phase flows. 

Once the testing procedures and equipment requirements were decided upon, 
efforts were made by both NASA and PFC to locate, obtain, and setup this equipment. 
PFC transformed their Mobile Clean Room into a supersonic nozzle testing facility to 
support these parametric studies. PFC fabricated and assembled a regulated supply 
system including regulator panels, water supply system, supply hoses, nozzle support 
stands, etc. B-CC supplied additional necessary equipment including lasers and optical 
benches. As stated above, NASA KSC supplied the nozzle and water injector assembly. 
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Measurements were taken over several days to obtain data describing the velocity 
profile of die nozzle using Breathing Air only. This was accomplished through die use of 
a hot-film probe, laser, and optical bench. The hot- film probe worked quite well, 
providing data that was consistent, repeatable and as theoretically predicted Using the 
calibration curve supplied by the probe manufacturer, die actual exit velocities could be 
calculated. The laser was used to maintain axial alignment of the probe with the nozzle. 
Measurements were taken over a range of pressures from 200 psig to 500 psig. Readi ng s 
were taken beginning at the exit plane of the nozzle and incrementally increased out to a 
distance of 10 inches away (axially) from the exit plane. 

During this phase of testing, observations of die flow characteristics of the nozzle 
clearly showed shock waves exiting the nozzle. Several black and white photographs 
were taken of these shock waves. The data gathered was analyzed and incorporated into 
the mathematical model as applicable. 

The next phase of testing involved introducing water droplets into the gas stream 
This was accomplished through the use of the NASA KSC designed injector assembly. 
This injector regulated die amount of water introduced into the system by of 3 
0 .W 25 inch orifice. System differential pressure of 20 psig forced the water through the 
orifice and into the gas stream. During initial setup and configuration of this gas-liquid 
p hase of testing die hot-film probe was inadvertently destroyed. The probe was 
determined to be irreparable and efforts were undertaken to devise alternative methods of 
data acquisition. 


Several different data acquisition methods were discussed between all parlies 
involved. The use of a Laser Doppler Velocimeter (LDV) was brought up as a viable 
approach to obtaining data on the gas-liquid system. The LDV has the capability to 
provide the desired information but it is an expensive piece of equipment that is not easily 
transported. High-speed photography (1,000,000 frames per second) was presented as a 
possible method of obtaining particle size and particle distribution data. Attempts have 
already been made by NASA to use such high-speed photography in the past, however 
wi the optics package that was available at file time, satisfactory results could not be 
oDtaincd. In addition, the photographic equipment and facilities could not be obtained 
wifinn the short time frame of this project. Possibly with a different optics package, this 
will prove to be a viable method of data acquisition in the future. 

It was determined that a pitot tube probe could be used to measure the desired 
velomty data. The pitot tube has the ability to withstand the detrimental effects of the 
liquid droplets in the system. The only drawback to the pitot tube probe is that it is 
somewhat intrusive to the flow being measured. 
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Efforts were made to locate and obtain a supersonic pitot tube probe. These 
efforts were unsuccessful. Therefore, with joint effort from both B-CC and PFC a 
suitable pitot tube probe was fabricated. The probe was a blunt end probe fabricated 
from a hypodermic needle and one-eighth inch copper tubing. This probe was coupled to 
a pressure transducer with strip chart recorder provided by NASA. 

The pitot tube probe was successfully used to measure velocities axially from the 
exit plane to 10 inches away. First die probe was used to obtain velocity data using only 
Breathing Air (no water). This data corresponded with the “gas only*’ data gathered using 
the hot-film probe earlier. Next the water droplets were introduced into die system as 
described above and testing continued. Again, measurements were taken over a range of 
pressures from 200 psig to 500 psig and from 0 to 10 inches away from the exit plane of 
die nozzle. The data gathered using the pitot tube probe was quite repeatable and as 
expected 

The flow characteristics of the nozzle with the gas-liquid mixture were observed 
during this second phase of testing. As before with gas only, shock waves were visible 

exiting the nozzle. These shock waves were again photographed using black and white 
film. 


All the data obtained from testing performed at PFC was analyzed by B-CC and 
ERAU. This data was incorporated into the mathematical model of die nozzle developed 
by ERAU. The model was “fine tuned** using this empirical data to insure that it did 
indeed correctly represent the existing nozzle configuration. 

Arrangements were made through LARC to use Laser Doppler Anemometry as 
another method to verify die two-phase flow data obtained with the pitot tube probe. The 
equipment required to perform such testing was located at NASA’s Langley Research 
Center. The Laser Doppler Anemometer (LDA) would be able to measure velocities of 
the gas-liquid stream exiting the nozzle without being intrusive to die flow. It was also 
thought that the LDA would be able to provide information on die size of the water 
droplets and well as die distribution of these water droplets. 

Members from KSC and B-CC traveled to LARC to discuss this third ph ag e of 
testing using die LDA. After the testing procedures and requirements were a gr eed upon, 
plans were made to conduct the testing. Members from both KSC and B-CC traveled to 
LARC to perform the testing. 
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The Laser Doppler Anemometer provided information on the velocity of the gas- 
liqmd stream at the desired locations along the axis of the nozzle. Team members were 
unable to obtain information regarding die size of the water droplets or the distribution of 
these water droplets. Even so the personnel on-site at LARC were of great assistance to 
us throughout this testing, often at great personal sacrifice to their schedules The 
velocity data obtained at LARC was analyzed by B-CC and ERAU and incorporated into 
the computer model as applicable. 


With the completion of the testing performed at LARC, the first objective of this 
project was met There now existed empirical verification of the theoretical design. The 
next step of the project was to determine a theoretical optimal nozzle design using the 
mathematical model. Unfortunately, acquiring testing equipment for use at PFC and 

D °PPler Anemometry equipment on-site proved to be 
logisticalfy difficult This led to several schedule delays and due to the short window of 
tne project made it necessary to revise the project requirements. 

The revision of die project requirements was discussed at the last status meeting 
Smce the tune allocated for the project was used up, it was decided that several items 
would be deleted from the Statement of Work. Items deleted include the requirement to 
mime the nozzle to its optimal configuration. Associated with this deletion is the 
requirement to fabricate and test any such newly designed nozzle configuration. 

i 0 » ** desired to use die computer model to determine several possible 

nozzle configurations that were dependent upon the specific application of the nozzle 
information would be coupled with a market study and commercialization plan with 
h^esto eventually market nozzle configurations/packages that proved feasible. These 
requirements were also deleted from this Statement of Work. 


A. th0U8 ?, Sevend req“ire»“Bts had to be deleted from die Statement of Work. 

ipr je« was still a success. The objective of obtaining empirical verification of the 
ftwened design was met The data showed that the water droplets did indeed exit with 
rop wsomet reloads. The development of a mathematical computer model representing 
*emasnng nozzle system was also accomplished. Although outside die scope of this 

nozzle ?“* con J uter modeI c “ he utilized to determine several different 

^^gMnnetiy s once die system parameters pertaining to each specific application arc 
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Once tiie computer model is available from ERAU, PFC plans to use it to 
investigate the feasibility of commercializing this nozzle technology. Initially, PFC 
hoped this characterization study would determine two things. First, what is the lowest 
supply pressure that will still provide supersonic exit velocities using the current n ozzl e 
design? Second, what is the smallest geometric configuration that will provide 
satisfactory cleaning ability? These questions were never answered during this study and 
therefore PFC will pursue these issues privately. The supply pressure issue is of concern 
to PFC because we feel that the required 300 psig could pose safety concerns if the 
nopde were to be used in commercial applications. Also, greater design flexibility can be 
enjoyed with smaller overall nozzle dimensions. 

If tiie computer model shows that it is theoretically possible to achieve the desired 
performance from these different nozzle configurations, then PFC will look into the 
fabrication and empirical verification of such nozzles. So far this study has been limited 
to the testing of a mixture of Breathing Air and demineralized water. Theoretical and 
empirical studies encompassing a variety of gas/liquid mixtures and system pressures 
would need to be conducted to verify nozzle performance once potential applications 
utilizing these different mixtures were identified. 

Once an application is shown to be technically feasible, a market study would be 
conducted to ascertain the need for such a product. If the market showed promise, a 
marketable “product” pertaining to that particular application would be developed. This 
“product” might range from developing and supplying a specific nozzle design per 
customer requirements, to developing, fabricating, assembling, and packaging a complete, 
self-contained, ready-to-use “nozzle system.” The nozzle could also be incorporated as 
part of a larger cleaning process. 


Initial applications of interest to PFC involve tile precision cleaning of K-bottle 
and similar cylinders, as well as utilizing the nozzle’s cleaning abilities in the electronics 
industry. PFC understands that a NASA designed K-bottle cleaner, incorporating the 
supersonic nozzle, already exists but as of yet PFC has not been able to perform any 
testing using this device as it has not been made available. PFC has already begun to 
compile a database of potential customers interested in using this type of nozzle 
technology in various commercial applications. 


PFC believes that several publications should come out of this characte rizati on 
study. The literature surveys performed by B-CC and ERAU could be published in 
literature survey journals. Articles written by B-CC and ERAU to specifically describe 
then* aspects of this study could be published in the appropriate technical journals. 
Articles written by PFC pertaining to the use of this nozzle technology in the cleaning 

industry could be published in the appropriate trade publications, such as “Precision 
Cleaning” magazine. 
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